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Duality violations in hadronic r decays and the value of a 
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Current estimates of the size of nonperturbative effects on the precision determination of 
a s from hadronic r decays are incomplete for at least three reasons: (i) In none of the existing 



C\| ' determinations has the effect of duality violations (DVs) been estimated quantitatively. We 
demonstrated in previous work [IJ |2] and in this talk that effects from DVs may turn out 
to be numerically significant, and that the use of only doubly-pinched moments of the r 
spectral functions in order to suppress DVs is not reliable, (ii) While the moments used by 
ALEPH[4] and OPAL[5] in principle probe the OPE up to dimension D=16, the analyses 
retain contributions only up to D=8. This was shown to be not self-consistent [3]. (iii) The 
most precise quoted values for a s are based on ALEPH data. However, it is now known that 
P-c correlations due to the unfolding of spectral functions were omitted in ALEPH's latest data 
o ,1 analysis [2]. Our preliminary explorations indicate that the impact of this on the error in a s 
is not necessarily negligible. 

The physical reason for the appearance of DVs is the presence of resonances in the vector 
(V) and axial (A) r spectral functions, which neither QCD perturbation theory nor the OPE 
describe adequately. As demonstrated in Ref . [1] , it is possible that the pattern of resonances 
is rather different in vector and axial channels, and this could lead to a significant effect on 
the value of a s even if the sum of V and A spectral functions (V + A) looks relatively flat 
in the region between s ~ 1 GeV 2 and s = m? T . In other words, the assumption that DVs 
can be safely neglected in V + A may turn out not to be justified. Because of this and point 
(ii) above, the estimates of the nonperturbative part of R T (and other moments) used in the 
literature need to be revisited. 

There is no quantitative theory of DVs in QCD, and in order to investigate this issue in 
more detail, one has to resort to a (physically reasonable) model [T]. This makes it possible 
to explore finite-energy sum rules with simple, nonpinched weights. Figure [T] shows a fit to 
the nonstrange integrated vector spectral function using the sum rule 

/•so ^ r f°° 

/ ds p exp (s) = -— (h ds TIope(s) - ds p D v(s) . (1) 
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Here p e xp{s) is the experimental spectral function, available up to sq = m^., Uope{s) the 
OPE expression for the vacuum polarization, and Pdv{ s ) the "duality- violating" part of 
the spectral function, i.e., the part not present in Im Uope{s + ie)- The contour integral 
over Uope{s) is parametrized by a s (m^.) and the OPE condensates. The functional form of 
Pdv( s ) is n °t known from first principles. Following Ref. p~|, we model it as 

Pdv(s) = k e~ ls sin (a + (5s) , (2) 
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Figure 1: Fit to Eq. (QP, with the ansatz (TJj] (red curve) and without the ansatz (green curve). 

which introduces four new parameters into the fits to experimental data. We restrict so to 
an interval Sq € [s m i n ,m^], and we vary s m i n in order to check for stability. For the physical 
motivation of the model of Eq. ([2]), see Ref. pp. The fit in Fig. [1] was done with OPAL data 
[5], and s min = 1.5 GeV 2 . The result for a s {m 2 ,) is 

a s (m 2 T ) = 0.307(18)(4)(5) (FOPT) , a s (m 2 T ) = 0.322(25)(7)(4) (CIPT) , (3) 

where errors are the x 2 error from the fit, from varying s m i n = 1.5 ± 0.1 GeV 2 , and from 
varying the estimated coefficient of (a s /ir) 5 in the Adler function in the range to 566. 

Fig. 1 demonstrates that DVs cannot, in general, be neglected in analyzing hadronic r 
decay data, even integrated versions thereof. This is confirmed by explorations which also 
include the axial channel in the fits. Essentially all existing r-based extractions of a s neglect 
DVs, and thus assume model (2) with k set to zero by hand. Such an assumption is necessar- 
ily accompanied by an uncertainty not accounted for in currently quoted errors for ct s . Fig. 
1 shows that it is dangerous to assume, a priori, that this uncertainty can be neglected for 
analyses involving doubly-pinched weights. We are presently studying different combinations 
of V and A spectral-function moments with the goal of quantifying this uncertainty. This 
may help reducing the fitting errors on a s , but at present, if one wishes to be conservative, 
and stick only to results not subject to additional, currently unquantified uncertainties, one 
should treat a s from r decays as being known to at best the accuracy given in Eq. ([3]). 
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